frep-ph/ 0006212 



June, 2000 
OU-HEP-352 



Polarized Structure Functions g2{x) in the Chiral Quark Soliton Model 

O 

O ■ 

O M. Wakamatsulil 

(N 

^ | Department of Physics, Faculty of Science, 



ON 



X 



Osaka University, 
Toyonaka, Osaka 560-0043, JAPAN 



> 

<N . PACS numbers : 12.39.Fe, 12.39.Ki, 12.38.Lg, 13.60.Hb, 14.20.Dh 

(N 

O ■ Abstract 
O 

The spin-dependent structure functions g2(x, Q 2 ) are investigated within the framework of 
the chiral quark soliton model. It turns out that the twist-3 part of g2(x, Q 2 ) gives nonnegli- 
gible contributions to the total distributions at the energy scale of Q 2 = 5 GeV 2 but mainly 
in the smaller x region only, so that the corresponding third moments Jq 1 x 2 g2{x, Q 2 ) dx are 



pretty small for both of the proton and neutron in conformity with the recent E155 data. 



In our opinion, the unexpectedly small quark spin fraction of the nucleon indicated by the 
EMC experiment |J and the light flavor sea-quark asymmetry revealed by the NMC measure- 
ment |§ are two remarkable discoveries in the field of nucleon structure function physics. Un- 
doubtedly, both are manifestation of nonperturbative QCD dynamics imbedded in the physics 
of high-energy deep-inelastic scatterings. An outstanding feature of the chiral quark soliton 
model (CQSM) is that it can explain both of these observations without recourse to any fine- 
tuning. (This potentiality of the CQSM was already noticed in @,f§].) In fact, we have already 
shown that it reproduces all the qualitatively noticeable features of the recent high-energy mea- 
surements, including the NMC data for F^x) — F 2 n (x), F^ix)/ F% (x) 0, the Hermes and E866 
data for d(x)-u{x) the EMC and SMC data for £(x), g?{x) and gf(x) ^0,0(13], in no 
need of adjustable parameters except for the starting energy scale of the renormalization-group 
evolution equation ||[6|,[/J. 
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Encouraged by this success, we now push on with our analyses to the twist-3 parton distri- 
butions in the nucleon. It is well known that, altogether, there are three twist-3 distribution 
functions - chiral odd, e(x, Q 2 ) and hi(x,Q 2 ), and chiral-even, g 2 (x,Q 2 ). They are gener- 
ally believed to provide us with valuable information on quark-gluon correlations in hadrons. 
In the present investigation, we shall focus our attention on g 2 (x,Q 2 ). There already exist 
several theoretical investigations of the polarized structure functions g2(x,Q 2 ). The most of 
those are based on various modifications of the MIT bag model as well as its original ver- 
sion [14], [L5|, [TBI, [T7|, [18| . There also exists an investigation based on the Nambu-Jona-Lasinio 



soliton model ||19|| , which is essentially equivalent to the CQSM. Confining to the lower mo- 



ments of g 2 (x, Q 2 ), also available are theoretical predictions based on the QCD sum rule lf20| , |21 



as well as the quenched lattice QCD |2"2" |, etc. 

Although based on essentially the same model as the present investigation goes far 
beyond the previous one in many respects. First, the polarization effects of the negative- 
energy Dirac-sea quarks in the hedgehog mean- field are fully taken into account. This is very 
important for offering any reliable predictions for antiquark distributions. Secondly, we also 
include the novel 1/N C correction (or the first order rotational correction in the collective 
angular velocity Q) to the isovector distribution functions. Without inclusion of it, some 
fundamental isovector observables like the nucleon isovector axial coupling constant would be 
largely underestimated, thereby being led to the so-called 11 g a. problem" in the hedgehog soliton 
model |2B|,|23|]. Thirdly, the nonlocality effects (in time) inherent in the theoretical definition 



of parton distributions are treated in a consistent way |5|,p5|. 

We start with the following definition of the distribution functions : 
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which is just standard except for the isospin dependence explained below. Here the parts 
containing 1 and r 3 respectively give isoscalar (/ = 0) and isovector (/ = 1) combinations of the 
relevant distributions. They are normalized such that the corresponding quark distributions 
<7<r (x) with q = u,d are given as 

9 { T /d \x) = \ i9 { T =0 \x) ± ff? =1) (*) ] (0 < x < 1), (3) 

and similarly for gi(x). The distribution functions (1) and (2) are formally defined in the 
region — 1 < x < 1. The functions with negative x are to be interpreted as giving antiquark 
distributions gx\x) with q = u,d according to the rule : 

9t'*\x) = \ ig { T =0 \-x) ± 9$ =1 \-x) ] (0 < x < 1) , (4) 



and similarly for g\ [x] . The corresponding structure functions for the proton and the neutron 
at the model energy scale are then constructed as 



g^\x) 



(5) 



The twist-2 distribution functions g\(x) were already evaluated in [EL0]. The distribution 
functions g^~°\x) and (x) can be evaluated within the same theoretical framework. 

Skipping the detail, here we only recall the fact that the isoscalar and isovector parts have 
totally different dependences on the collective angular velocity Q of the rotating hedgehog 
mean field, which itself scales as l/N c : 



g^ =0 \x) ~ N c O(tt l ) ~ 0(N° C ), 

g¥=°\x) ~ N c [0(n°) + O^ 1 )] ~ 0(iV c 1 ) + O(N c ) . 
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Figure 1: (a) The isoscalar parts of the distribution functions gi(x) and gr(x). (b) The 
isovector parts of gi(x) and gr(x)- 

Shown in Fig.l are the results for gx~°\x) and g^~ (x) in comparison with the twist- 
2 distributions g[ J ~°\x) and g± (x). We first point out that the isoscalar and isovector 
distributions have totally dissimilar shapes reflecting quite different iV c -dependence given in 
(6) and (7). Next, comparing g^ (x) and g[ I ~°\x), one finds that the g^~°\x) has a peak at 
smaller value of x than g[ I ~°\x) and damps faster as x increases. The same tendency is also 
observed for the isovector distributions, but in this case the concentration of the distribution 
g? (x) into the smaller x region is even more profound. This is due to more significant 
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effect of vacuum polarization, which is peaked around x ~ 0. It means that the "valence- 
quark-only" approximation adopted in |TJ| cannot be justified at least for g!£~ (x) and to a 
lesser extent also for gi (x). Note, however, that the vacuum polarization contributions, 
which are dominant in the smaller x region, are partially canceled in the corresponding spin 
structure function 02 0*0 defined as a difference of g^~ 1 \x) and g[ I ~ 1 \x). In any case, 
we emphasize the following. The CQSM gives fairly different predictions for the shapes of 
g\{x) and gr{x)- Moreover, the shapes of both distributions are strongly dependent on the 
isospin (or more generally flavor) combinations. Furthermore, both distributions gi~ l \x) and 
^(x) have large support in the negative x region, which implies sizable flavor asymmetry 



of the spin-dependent sea-quark (antiquark) distributions jq,|26j. 

To compare these predictions of the CQSM with the existing high-energy data for g%{x, Q 2 ) 
and g$(x, Q 2 ), we must take account of the scale dependence of the distribution functions. We 
have done it as follows. Remember first that the spin structure function #2 (a?, Q 2 ) is defined as 
a difference of gx(x, Q 2 ) and gi(x, Q 2 ), i.e. g 2 (x,Q 2 ) = gT{x,Q 2 ) — gi(x, Q 2 ). (In the following 
discussion on the scale dependence, we omit the isospin indices for the structure functions, for 
notational simplicity.) The Burkhardt-Cottingham sum rule p7] holds exactly, i.e. 



/ g2{x 1 Q 2 )dx = 
Jo 



if the charges (or the first moments) of gi(x, Q 2 ) and gr{x, Q 2 ) are equal, which in turn follows 
from the rotational invariance of the whole theoretical scheme. This property is automatically 
satisfied in usual low energy models like the MIT bag model 0,|IJ,[T§,[T7|,|18| or the CQSM |IJ . 



The g2^x, Q 2 ) is further decomposed into the twist-2 (Wandzura-Wilczek) part and the genuine 
twist-3 part as 



i/2 



x,Q 2 ) = - 9l (x,Q 2 ) + f 1 ^- 9l ( yi Q 2 ) (9) 

Jx y 

H^Q 2 ) = 92(x,Q 2 ) - gf w {x,Q 2 ). (10) 



For the QCD evolution of the structure function gi(x, Q 2 ) and twist-2 piece of giix-, Q 2 ), the 
ordinary Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation can be used. Here, we 
use the leading-order Fortran code provided by Saga group J2^J. The flavor non-singlet and 
singlet channels are treated separately, since the mixing with the gluon distribution occurs 
in the latter. On the other hand, the Q 2 -evolution of twist-3 distributions is known to be 
quite complicated due to mixing with quark-antiquark-gluon operators, the number of which 
increases rapidly with spin or the moment of the distributions. However, Ali, Braun and Hiller 
found that, in the large N c limit, the Q 2 -evolution of the chiral-even twist-3 flavor-nonsinglet 
distribution gi{x-, Q 2 ) is described by simple DGLAP type equation with slightly different forms 
for the anomalous dimensions from the twist-2 distributions [3(J. Accordingly, the moments 
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of g2(x, Q 2 ) obey the following simple equation : 



M n [g 2 (Q 2 )} = L^M n [UQ 
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where M n [g{Q 2 )} = J* dx x™" 1 g(x, Q 2 ), L = a s {Q 2 ) / 'a s (Q 2 ni ), b = f N c - § N f , and 



:i2) 



with S n = YJj=\ ^- The Q 2 -evolution of the corresponding distribution functions can be 
handled by the method described in pifl . In principle, the flavor-singlet part of g~2(x,Q 2 ) 
mixes with the gluon distribution, and the Q 2 -evolution of it is not given by a simple equation 
as above even in the large N c limit. In the following study, we shall neglect this mixing 
effect with gluons in the twist-3 flavor-singlet distributions, for simplicity. Finally, the total 
g 2 {x-, Q 2 ) at the desired energy scale is obtained after combining the twist-2 and twist-3 pieces 
of g2(x,Q 2 ), which are evolved separately. 
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Figure 2: The theoretical structure functions x 2 g^x.Q 2 ) and x 2 g$(x, Q 2 ) at Q 2 = 5 GeV 2 
(solid curves) are compared with the corresponding E143 and E155 data (respectively shown 
by crosses and filled circles) as well as the predictions of Song's center-of-mass MIT bag model 
(dashed curves). The twist-2 parts of the theoretical structure functions (dash-dotted curves) 
are also shown for comparison. 



We show in Fig. 2 the theoretical structure functions g2(x, Q 2 ) at Q 2 = 5 GeV 2 for the proton 
(a) and the deuteron (b) in comparison with the corresponding experimental data. Here, the 
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crosses and the filled circles respectively stand for the E143 [32| and E155 data p3 |. The final 
theoretical predictions for g 2 (x,Q 2 ) are represented by the solid curves, while their twist- 2 
part g 2 {x, Q 2 ) are shown by the dash-dotted ones. The predictions of the center-of-mass 
MIT bag model by Song are also shown for the sake of comparison |[L8| . We point out that 
the predictions of the CQSM for both of g 2 (x, Q 2 ) and g 2 { x i Q 2 ) are relatively close to those 
of another version of MIT bag model given by Strattmann shown in j32| and sizably larger 
than those of Song's results. We also find that, according to the predictions of the CQSM, the 
differences between the full g 2 (x, Q 2 ) and their twist-2 parts are relatively small except for the 
smaller x region, although it is not clear from Fig.2 in which x 2 times g2{x,Q 2 ) are plotted. 
This tendency is also close to Stratmann's results rather than Song's results. 
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Figure 3: The twist-3 parts of the structure functions, x 2 g 2 (x, Q 2 ) and x 2 g 2 (x,Q 2 ) are com- 
pared with Song's predictions as well as the corresponding experimental data from his paper. 



To see it more clearly, we show in Fig.3 our results for the twist-3 part of g 2 (x, Q 2 ) and 
g 2 (x, Q 2 ) in comparison with Song's predictions. Here, the solid curves are the predictions 
of the CQSM, while the dashed curves are those of the center-of-mass bag model by Song. 



The experimental data in this figure are from ||18|| . One clearly sees that the twist-3 parts of 
g 2 (x,Q 2 ) are much smaller in the CQSM than in Song's bag model calculation. Because of 
the large uncertainties of the available experimental data, it is difficult to say at the present 
moment which theoretical prediction is favored. Nonetheless, small twist-3 contributions to 
the spin structure functions g 2 (x, Q 2 ) appears to be favored by the recent E155 analysis of the 
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twist-3 matrix element : 
d 2 {Q 2 ) 



3 / x 2 g 2 (x,Q 2 )dx 
Jo 



x 2 [gi(x,Q 2 



dx . 



(13) 



(Also interesting to notice here would be the fact that the instanton-liquid model of the QCD 
vacuum offers a qualitative explanation of the suppression of the twist-3 matrix element of 
g 2 (x,Q 2 ) relative to the twist-2 one We compare in Fig.4 the predictions of various 

theoretical calculations with the recent E155 data. As already shown in ||33|| , the predictions 
of some models are apparently incompatible with the E155 analysis, although we must be 
cautious about difficulties in obtaining reliable and precise experimental information for these 
quantities. For instance, large and negative d 2 for the neutron predicted by the QCD sum 



rules |20],|21[ apparently contradicts the E155 data. Similarly, large and negative d 2 for the 
proton predicted by the lattice QCD p2f seems incompatible with the E155 data. 
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Figure 4: The predictions of various theoretical models for the twist-3 matrix element d 2 
for the proton (a) and the neutron (b) are compared with the recent E155 analysis. Shown 
theoretical models are from left to right : QCD sum rules [20|, 2l|, lattice QCD |fffi| , MIT bag 



models MM, and the CQSM. 



On the other hand, the predictions of the CQSM as well as those of the MIT bag models 
seems consistent with the E155 data at least qualitatively. We see that the prediction of the 
naive MIT bag model for d\ is accidentally close to that of the CQSM. Its prediction d 1 ^ = 
also seems to lie within the experimental error bars. Note however that the predictions of 
the naive MIT bag model must be taken with care, since the SU(6) structure of the bag wave 
function (this is the cause of the result d^ = jlfj) apparently contradicts large and negative 
behavior of the twist-2 neutron structure functions gi(x,Q 2 ) confirmed in several previous 
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experiments [T^,[TT|,|T^,[T3|. This shortcoming of the original MIT bag model is remedied in 
Song's modified one, in which sizable SU(6) symmetry breaking effects are incorporated by 
hand, so that it reproduces both of gf(x, Q 2 ) and g™(x, Q 2 ) JIB| . However, now we sees that the 
MIT bag model so refined gives a prediction for d?,, which is about two times larger than that 
of the CQSM and lies outside the errorbars of the E155 data. More precise experimental data 
are absolutely awaited for drawing more decisive conclusion about the twist-3 contributions 
to the nucleon spin structure functions, thereby selecting various models of nucleon internal 
structure. 

To sum up, it has been shown in a series of paper that the CQSM reproduces all 

the qualitatively noticeable features of the recent high-energy data for the twist-2 structure 
functions of the proton, the neutron and the deuteron, with no adjustable parameter except 
for the ititial energy scale of the DGLAP evolution equation. In the present investigation, we 
have extended this parameter- free analyses to the twist-3 spin structure function g 2 (x, Q 2 ) . The 
theoretical predictions are shown to be consistent with the E143 and E155 measurements for 
g\{x, Q 2 ) and g$(x, Q 2 ) at Q 2 = 5 GeV, although the uncertainties of the existing experimental 
data are still too large to draw a decisive conclusion. We have also shown that the CQSM 
predicts very small twist-3 matrix elements d 2 for the proton and the neutron in conformity with 
the recent E155 analysis. The accumulation of more precise experimental data for g%{x, Q 2 ) as 
well as gi(x,Q 2 ) is absolutely necessary for more complete understanding of the nucleon spin 
structure. 
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